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Resumo 

Neste trabalho, recorremos a metodologia da transformada de Laplace a fim de mostrar a 
sua importancia na abordagem de uma classe de equagoes diferenciais fracionarias. Em par¬ 
ticular, apresentamos aplicagoes desta metodologia ao discutirmos possiveis generalizagoes de 
certos problemas fi'sicos no campo da viscoelasticidade linear e osciladores harmonicos, com- 
provando que o uso do calculo fracionario em modelagem e resolugoo de problemas usualmente 
abordados pelo calculo de ordem inteira, oferece vantagens promissoras para nos fornecer for- 
mulagoes mais consistentes com os dados experimentais. 

Palavras-chave: Calculo fracionario, Equagoes diferenciais fracionarias; Transformada de La¬ 
place. 

In this paper, we resort to the Laplace transform method in order to show its efficiency 
when approaching some types of fractional differential equations. In particular, we present 
some applications of such methods when applied to possible generalizations of certain physi¬ 
cal problems in linear viscoelasticity and harmonic oscillators, proving that fractional calculus 
is well suited for the modelling and solving of problems usually treated by ordinary integer 
calculus, with the promissing advantages of being able to provide more accurate theoretical 
predictions to fit with experimental data. 

Keywords: Fractional calculus, fractional differential equations, Laplace transform. 


1 Introdugao 

O calculo de ordem nao inteira, popularmente conhecido como calculo fracional ou Calculo Fra¬ 
cionario (CF), de uma maneira simples, tern a intengao de generalizar o calculo integral e diferen- 
cial, conforme proposto, independentemente, por Newton e Leibniz. Aqui, evitamos, sempre que 
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possfvel, o aparato matematico envolvendo explicitamente as formulas advindas, por exemplo, das 
diversas maneiras de calcular uma derivada. Ainda mais, utiliza-se a nomenclatura CF por entender 
que o nome esta totalmente consolidado e, em lingua portuguesa, e uma tradugao livre de fractional 
calculus. Um estudo versando sobre a linha do tempo envolvendo o CF pode ser encontrada em 
[mi enquanto um capi'tulo sobre a historia do calculo fracionario pode ser encontrada em 

Como ja mencionado, existe mais de uma maneira de calcular a derivada e, portanto, parece emi- 
nente questoes do tipo: Para que serve o CF? Onde utiliza-lo? Existe uma interpretagao geometrica 
e/ou fisica? Qual a relagao, se e que existe, com o calculo de ordem inteira? Muitas outras questoes 
de ordem matematica ou de ordem fisica podem ser colocadas. Por exemplo, apenas para men- 
cionar duas dessas questoes: Sera que existe uma regra da cadeia associada a derivada de ordem 
fracionaria? E o teorema fundamental do calculo, teorema que coroa os calculos diferencial e inte¬ 
gral, tern um analogo fracionario? Ainda mais, no calculo de ordem inteira emerge uma classe de 
fungoes, as chamadas fungoes especiais, solugoes das equagoes diferenciais ordinarias e/ou parciais 
que descrevem um particular sistema, enquanto no CE emerge tambem uma classe de fungoes a 
ele associado, solugao de um particular problema cuja derivada e de ordem nao inteira, as fungoes 
especiais do CE EiiziiHiiiiiniiiiiiiiiiiiiiiiiKiiiiiziiiiiiiiin]- 

O objetivo principal desse trabalho e introduzir as ferramentas basicas a fim de discutir um 
particular problema, composto por equagao diferencial fracionaria e condigao inicial. De modo a 
atingir este objetivo, isto e, discussao e resolugao de uma equagao diferencial fracionaria, aborda- 
se o tema derivada fracionaria que, por sua vez, requer o conceito de integral fracionaria. De 
uma maneira bastante simplificada, introduzimos o conceito de integral fracionaria para depois, 
utilizando tal conceito, introduzir a derivada fracionaria que, nesse trabalho discute apenas as 
formulagoes conforme propostas por Riemann-Liouville e por Caputo 01 inni na [H na [m |22] . 
Em particular, apesar de mais restritiva, a derivada de Caputo admite a mesma interpretagao para 
as condigoes iniciais que a formulagao classica de ordem inteira [III na ng. Mostrou-se tambem 
que essas duas formulagoes podem ser recuperadas a partir da definigao de Griinwald-Letnikov [20] . 

De uma maneira simples e objetiva, pode-se pensar nos operadores de ordem fracionaria como 
os operadores que representam fungoes da memoria sobre a historia de alguns sinais de sistemas 
ffsicos. Por exemplo, uma integral de primeira ordem, de uma variavel que representa o estado de 
um sistema, pode ser pensada como uma soma em que se atribuem pesos a cada ponto, com todos 
os pontos ponderados com o mesmo peso, independentemente de quao longe estao no passado. 
Isso sobre toda a historia do sistema. Uma integral de ordem fracionaria e tambem uma soma 
ponderada, mas com os pesos diminuindo para tras no tempo [HimiisiiTi]. 

Em termos matematicos, o CF atrai por si proprio um grande interesse, pois o formalismo en- 
volvido faz uso de diversas fungoes especiais, tais como a generalizagao da fungao fatorial, que e 
a fungao gama e a fungao de Mittag-Leffler de um parametro, como uma generalizagao da fungao 
exponencial, dentre muitas outras de interesse academico e pratico Existe mais de uma 

formulagao possfvel para o CF, sendo cada uma dessas mais adequada a um certo contexto ffsico do 
que outro. As definigoes mais comuns sao as de Riemann-Liouville, Caputo, Criinwald-Letnikov, 
Liouville, Weyl e Riesz-Feller [20] . Parece que o mimero de definigoes nao para de crescer [21] . Sao 
recentes as formulagoes de Hilfer [7], para particulares valores do parametro associado a ordem da 
derivada, onde as formulagoes de Riemann-Liouville e Caputo sao casos extremes e a formulagao 
proposta por Khalil-Horani-Yousef-Sababheh [23] a assim chamada conformable fractional deriva¬ 
tive que, numa conveniente tradugao para o portugues, pode ser chamada de derivada fracionaria 
compatfvel. 

Neste trabalho discutem-se modelos fracionarios associados a derivada de Riemann-Liouville e 
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de Caputo, no sentido de que tais modelos sao mais representatives que os respectivos modelos 
lineares classicos de ordem inteira, isto e, descrevem com maior acuracia os sistemas em questao. 

O trabalho esta disposto da seguinte forma: Na segao dois sao introduzidos os operadores de 
integragao e os operadores de diferenciagao fracionarios conforme propostos por Riemann-Liouville 
e Caputo. Na terceira segao sao abordadas as equagoes diferenciais fracionarias (EDF) conforme 
as formulagoes supra-citadas, sendo que para estas duas formulagoes apresentam-se, atraves de 
teoremas, os cases gerais. Na segao quatro, discute-se a transformada de Laplace associada as 
integrals e derivadas fracionarias e justifica-se a conveniencia de se trabalhar com a formulagao da 
derivada segundo Caputo. Ainda na segao quatro apresenta-se, atraves de exemplos, o calculo da 
transformada de Laplace de fungoes de Mittag-Lefller, concluindo a segao, resolve-se duas EDFs 
atraves da transformada de Laplace, cuja solugao e dada em termos de fungoes de Mittag-Leffler de 
dois parametros. Na segao cinco discute-se a modelagem de dois sistemas onde o calculo fracionario 
desempenha papel preponderante. Apos uma revisao do conceito de viscoelasticidade linear, em 
particular, discutindo os classicos modelos de Maxwell e Voigt, apresenta-se o modelo fracionario de 
Scott-Blair o qual, nos limites extremes, recupera os dois modelos classicos. O modelo fracionario 
tern sua justificativa atraves de graficos a partir dos quais fica clara a importancia da formulagao. 
Por fim, o classico problema do oscilador harmonico em sua versao fracionaria e discutido sendo 
a solugao dada em termos de fungoes de Mittag-Leffler com um e dois parame-tros m- Aqui 
tambem, o caso limite (classico) e recuperado no caso em que o parametro associado a derivada e 
igual a dois. 


2 Operadores de Integragao e Diferenciagao 

Visto que a definigao de derivada de ordem nao inteira nas formulagoes de Riemann-Liouville e 
Caputo dependem da integral fracionaria, comegamos com tal conceito, isto e, introduzimos a 
integral fracionaria para depois apresentar a derivada fracionaria. 

Dentre as diversas formulagoes existentes para os operadores de integragao e diferencia-gao 
fracionarias [131211 HU, abordamos, para os propositos deste trabalho, apenas as versoes segundo 
Riemann-Liouville e segundo Caputo, definidas a seguir. 

Definition 1 Sejam = [a,6] C K um intervalo finito e f € Li [a, 6]. As expressoes 6 ^b-f 

definidas por: 

(Xa+Z) (a^) = f {t)dt, (1) 

com X > a, v > 0 e 

(It-f) (x) = {t - f (t) dt, (2) 

com X < b, u > 0, onde T {n) e a fungao gama, definem as integrais fracionarias de Riemann- 
Liouvill^ (IFRL) de ordem n G M. num intervalo real finito. As integrais em Eq. m e Eq.([3) sao 
chamadas de integrais fracionarias a esguerda e a direita, respectivamente. 


^Alguns autores ElITo], distinguem a nomenclatura, baseandose nos limites inferior e superior das integrais. 
Segundo suas nomenclaturas, Eq. {TJ e Eq.© sao chamadas de versao segundo Riemann; e quando a = —oo e b = oo, 
denotam por versao segundo Liouville. 
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Definition 2 Sejam = [a, &] C K. um intervalo finito e f G AC’^ As expressoes e 

T^b-f definidas par 


com X > a e 


{v:+f) {x) = [(!"+"/) (a:)] 

r /wdf 

r(n-z^)A 


{T^b-f) i^) 


-Da^ [{ir''f) (^)] 

(-^)" f jt) 

r(n-z^)7, (t-:r)^+" 


( 3 ) 


( 4 ) 


com x < b, onde n = [i/] + 1 e [v] e a parte inteira de v, definem as derivadas fraciondrias de 
Riemann-Liouville (DFRL) de ordem ^ S K, com p > 0, d esquerda e d direita, respectivamente. 


Definition 3 Sejam = [a, 6] C K e G M, com p > 0. Considere ojDx ^ x^b o,s DFRL como 
nas Eq. e Eq.(|4]) e defina n = [i/] + 1, ^ No; n = z/ se z/ S No, entao as expressoes 


C-r^u 


V’jfix) =a VI 


/c=0 


-(a: — of 


( 5 ) 


e 


^,VlJ(x) =x VI 




definem as derivadas fraciondrias de Caputo d esquerda e a direita, respectivamente. 


( 6 ) 


Observamos inicialmente, que todas as definigoes acima se reduzem aos casos classicos quando 
escolhemos um valor inteiro para a ordem o. Alem disso, as definigoes que apresentamos foram 
elaboradas para um zz G M, no entanto, elas sao igualmente validas se considerarmos p G C, com 
Re(z^) > 0, [HllIH]. 


3 Equagoes Diferenciais Pracionarias 

Nesta segao apresentamos as chamadas EDF e alguns resultados sobre a existencia e unicidade 
das solugoes, quando definidas, em intervales finitos. Observamos que existem varies estudos com 
resultados eventualmente distintos para o problema de existencia e unicidade de solugoes destas, 
isto porque cada formulagao depende das definigoes usadas para os operadores de diferenciagao e 
integragao 0 [13 [H [ig. Os problemas mais estudados envolvem as solugoes de EDFs segundo 
Riemann-Liouville e segundo Caputo: o primeiro por ser a definigao mais difundida e o segundo 
pela conhecida interpretagao fisica das condigoes iniciais (ou de fronteira). 
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3.1 Formulagao Segundo Riemann-Liouville 

Uma equagao diferencial (nao linear) fracionaria de ordem z/ > 0, definida num intervale finite [a, b] 
e da ferma 

{V';^+y)ix) = fix,y{x)), (7) 

onde e o eperader de DFRL. 

Se estabelecermos um cenjunto de cendigoes iniciais (eu de frenteira) 

{K^^-^y) («+) = lim (K-'^-^y) (x) (8) 

x—¥a'^ 

= fc = 0,1,... ,n - 1, 

onde € M, n = [z/] + 1 se z/ ^ N e z/ = n se z/ € N, entao analogamente ao case de ordem 
inteir£0, denotamos a Eq.d?]) junto com as condigoes Eq.® de um problema de Cauchy (fracionario). 
Chamamos atengao para o caso em que fc = n — 1 nas condigoes acima, pois devemos interpretar 

{v:--y) (a+) = lim (P^-y) (x) 

x—¥a'^ 

= lim (X,V 2 /)(^)- 

Foge do escopo deste trabalho a discussao da existencia e unicidade da solugao do problema 
de Cauchy fracionario, entretanto, mencionamos que tais resultados ja se encontram formulados 
na literatura e, podem ser encontrados em 01111121], por exemplo. Ressaltamos tambem que o 
problema de Cauchy (vide Eq.([7]) e Eq.®) pode ser formulado, equivalentemente, em termos de 
uma equagao integral de Volterra Emm]. 

Para os propositos deste trabalho, entretanto, a formulagao do problema de Cauchy como ex- 
plicitado acima e muito geral, assim nos restringimos a discutir exemplos para o caso em que as 
respectivas equagoes diferenciais ordinarias lineares fracionarias (EDOLF) sao da forma 

V^+yit) - Xy(t) = f(t) (9) 

com A S M e as condigoes iniciais 

W] fc = 0,1,..., n - 1, (10) 

onde n = [z/] + 1. Observamos ainda que o limite inicial de integragao foi escolhido como sendo o 
ponto t = 0, pois mesmo que tenhamos inicialmente um outro ponto t = a, sempre e possivel por 
meio de uma mudanga de variavel fazer uma translagao para a origem. 

A solugao deste problema que, em geral, e conduzido a uma integral de Volterra apresenta 
solugao em termos de uma integral envolvendo uma fungao de Mittag-Leffler com dois parame-tros 
Ea,p (Oi ^ explicitamos a seguir, para referenda futura. 


yit) = Y, (xn 

(A (t - r)'") / (r) 


k^O 

pt 


dr. 


(11) 


dando-nos uma solugao explicita para a equagao integral de Volterra associada ao problema de 
Cauchy composto pelas Eq. (0) e Eq. (ITUl) 0 [HI HI] . 


^Obtido quando i/ G N. 
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3.2 Formulagao Segundo Caputo 

Analogamente a formulagao segundo Riemann-Liouville, consideramos uma equagao diferencial (nao 
linear) fracionaria de ordem i/ > 0, definida num intervalo finito [a, b] da forma 

iaT^xV) (x) = f ix,yix)), (12) 

onde a'^x 6 o operador de DFC, sujeita as condigoes iniciais 

(fp^2/)(0) = 6fe, fc = 0,l,...,n-l, (13) 


onde n=[i/] + lse^^Ne^ = 
segundo Caputo, as derivadas 


n se ^ G N. Observamos que na formulagao do problema de Cauchy 


C 

a 


V 


k 

X 


dfc 

da;^ 


sao as proprias derivadas usuais de ordem inteira, o que nos leva a tradicional interpretagao fisica 
das condigoes iniciais, como no usual problema de Cauchy para ordens inteiras. 

O teorema de existencia para o problema de Cauchy pode ser encontrado em [22] e, equiva- 
lentemente a formulagao de Riemann-Liouville, este problema pode ser conduzido a uma integral 
de Volterra (nao-linear) [501 HI] • 

Nessa formulagao, o teorema de existencia para o problema de Cauchy supracitado fica entao 
enunciado da seguinte forma [22] . 

Ressalta-se que os resultados entre as formulagoes de Riemann-Liouville e de Caputo, sao tais 
que enquanto na primeira garante uma solugao continua somente em (0, h] na segunda a solugao e 
garantida ser continua em [0, h] [501155] . 

Novamente, para os propositos deste trabalho as formulagoes acima sao muito gerais, assim nos 
restringimos a discutir exemplos para o caso das EDOLF da forma 


oKy(.t) - = fit), 


(14) 


com A G R, junto com as condigoes iniciais 

[oT^xVit)] (^0 = fc = 0,1,..., n - 1. (15) 

E, novamente, uma vez conduzida a sua forma integral, e possivel obter uma solugao explicita 
|H] para o problema de Cauchy composto pelas equagoes Eq.(|Ill) e Eq.dH]) em termos das fungoes 
de Mittag-Leffler como 


yit) = ^ bkt''E„^k+i (At'") -f 




(A (t - tY) f (t) 


(t-r) 


1 - 1 ^ 


dr. 


(16) 


4 Transformada de Laplace dos Operadores Pracionarios 

Sabemos do calculo de ordem inteira que o metodo da transformada de Laplace e uma ferramenta 
muito litil na analise de equagoes dife-renciais lineares, principalmente no caso em que a equagao 
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possui coeficientes constantes. Neste caso, a transformada de Laplace reduz a EDO numa equagao 
algebrica que e, em geral, muito mais simples de se solucionar, deixando a dificuldade de se obter 
a solugao final da EDO de partida a um problema de inversao. Nesta segao, veremos que a mesma 
metodologia tambem pode ser usada para resolver problemas de valor inicial relacionado com as 
EDOLF (vide Eq.® e Eg. lfTHl i. 

Lembremos que uma fungao / (t) e dita de ordem exponencial a se existir constantes positivas 
M eT, tais que 

e-“‘|/(t)| <M, 

para todo t > T. Dessa forma, dada uma / : [0, oo) —>■ M, de ordem exponencial a a fungao F 
definida por 

POO 

F{s) = / (17) 

Jo 

e chamada a transformada de Laplace da /, com a condigao Re (s) > a que garante a existencia 
da integral acima. Da nomenclatura classica e usual denotar a transformada de Laplace de / por 
C[f{t)] = F{s) e iremos denotar por [-F'(s)] = fit) a transformada de Laplace inversa, que 
como sabemos e unica pelo teorema de Lerch |25] . 

Dentre as diversas propriedades conhecidas da transformada de Laplace, listamos as seguintes: 








Linearidade: C [af{t) + (3git)] = aC [f{t)] + pC [git)] = aF(s) + pCis), com a,P gC; 

Convolugao: £[/(<)*g(<)] = C[fit)]C[git)] = F(s)G'(s), onde o produto * e o de con- 
volugao; 

Derivadas: 


d"/ 

dt" 


n — 1 

= s^Fis) - 


• Integral: 


£ 


[ f (t) dr 
Jo 



Assim, nosso intuito e calcular as transformadas de Laplace da IFRL, da DFRL e da DEC de 
uma dada fungao / (t) suficientemente bem comportadc@. 

Lembrando que a IFRL pode ser escrita como o produto de convolugao 


iK+f) it) = V y) 

= fit)* Pu it), 


entao segue da transformada do produto de convolugao listada acima, que 

£[(J„%/) it)]=C[fit)^p, (t)] 

= ns)\- 


( 18 ) 


^No sentido de satisfazer as condigoes de existencia das IFRL/DFRL e das respectivas transformadas. 
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Agora da Definigao[2j temosque (X>o+/) (0 = (t)^ com n— onde g{t) = (t), 

assim 


c [{v^o+f) ( 0 ] = ^ [K+git)] 


n-1 

= s"G(s) - 


(19) 


Mas, por outro lado, temos 


G(s)=£[(l"-‘'/) W] 

_ F{s) 

gti-j/ ’ 

portanto, usando a Ea. (|2ni) na Ea. lfT^ . obtemos 


G m+f) W] = s'^Fis) - E 




onde = lini^ 2?^ (Xq_|_ “"/) (i). 

Por outro lado, quando usamos a Definigao [3l podemos escrever que 

= ^G(s), 

s^-’' ^ ' 


( 20 ) 


( 21 ) 


( 22 ) 


onde 


G(s) = C [(Xo"+/) (t)] 

n-1 

= s"F(s) - E (0+). (23) 

Logo, os resultados das Eq. (P^ e Eq. (P5|) nos levam a seguinte expressao 

n-1 

F (t)] = s'^Fis) - E 5^'-'=-(0+). (24) 

fc=0 

Novamente fica evidente a diferenga entre as definigoes dos operadores de DERL e DEC, sendo 
que a segunda se mostra nitidamente mais oportuna de ser usada quando temos interpretagoes 
fisicas claras das condigoes iniciais do problema de Cauchy, como nos mostram os resultados das 
transformadas de Laplace para a DERL e DEC (vide Ea. d^ e Ea. (l24l) l. 

Antes de apresentarmos alguns exemplos, listamos, por conveniencia, as transformadas de La¬ 
place das fungoes de Mittag-Lefflei@, que nos serao uteis na hora de resolvermos o problema de 

®Estas transformadas podem ser facilmente calculadas, usando a representagao em serie das fungoes de Mittag- 
LefHer e calculando as transformadas termo a termo, o que e possivel visto que estas fungoes sao inteiras no piano 

C. 
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inversao C ^ : F{s) f{t). Assim, temos 


c[E^i-xn] = 


ga-1 

s“ + 1 
1 + As-“ 




c 




s“ + A 
1 + As-“ 

gOtJ-P 

(s“ + A)'^ 
s~^ 

(1 + xs-°^y 


(25) 


(26) 


Antes de passarmos a discutir modelos via CF, apresentamos dois problemas de 
envolvendo a formulagao de Riemann-Liouville (um problema homogeneo e o outro nao 
enquanto a formulagao de Caputo sera aprensentado na proxima segao. 

Example 4 Determine y(t) no problema 

I D§g_y{t) - Xy{t) = 0, 

\ p|+-'y(0+) = l|+2/(0+) = 6o G K. 

Aplicando a transformada de Laplace a equagdo, obtemos 


(27) 

valor inicial 
homogeneo) 


o que implica que em 


donde podemos verificar, usando a Eq. d^S)) . que 

y{t) = bof-^Ei i{Xt^), 


C ^D^+yit) - Aj/(t)J = 0 
s^Y{s) — bo — XY (s) = 0 
_ xj Y{s) = bo, 


YA = 

52 — A 


onde ^a ,/3 (*) e o, fungdo de Mittag-Leffter com dois pardmetros m- ^ solugdo coincide eom a 
solugdo do problema de Cauchy conforme a Eq. CD- 

Vejamos ainda um segundo exemplo, agora nao-homogeneo. 

Example 5 Resolva o problema composto pela EDF e as condigoes a seguir: 

TX^+yit) - Xy{t) = C, 

'dI-\{Q+) = 4+j/(0+) = 6o G R, 

I?|+\(0+)=4+y(0+) = 6iGM. 
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Aplicando a transformada de Laplace na equagdo, obtemos 


ou ainda 


0 que implica em 


- Xy{t) 
r(3) 


r(3) 


3 — y(s) — ——h sbi + 6o, 


Y{s) = 


sbi r(3) 1 


4 . ' 4 , ' Q 4 , ? 

S3— A S3— A s S3— A 


donde podemos verificar, usando a Ea. (l26l) e a propriedade da transformada do produto de con- 
volugao, que 


y{t) = bot^Ei i (^^0 


+ 

= y^bkt 


4-fc-ip 

3 ±!/4 4 

3 ’ 3 


/c=0 
rt ^ 


(«») 




^dr 


e, novamente, o resultado coincide com a formula para a solugdo do problema de Cauchy conforme 

a Eg. (ITT]). 


Estes exemplos apesar de simples, ilustram que a metodologia da transformada de Laplace con- 
tinua tendo a mesma eficacia para a resolugao de EDOLFs (com coeficientes constantes) para a 
formulagao de Riemann-Liouville, assim como nos casos de ordens inteiras. Conforme menciona- 
mos,uma aplicagao usando a formulagao segundo Caputo, sera discutida quando apresentarmos 
uma possivel generalizagao para o pro-blema de Cauchy associado aos osciladores harmonicos na 
Segao 15.21 


5 Modelos via Calculo Fracionario 

Vejamos agora como os conceitos apresentados do calculo fracionario podem ajudar a descre-ver 
de forma mais apropriada certos fenomenos fisicos, quando comparados com a abordagem classica. 
Ilustraremos a afirmagao acima investigando um problema particular no campo da viscoelasticidade 
linear e numa possivel genera-lizagao do problema associado ao oscilador harmonico. 

5.1 Viscoelasticidade Linear 

A viscoelasticidade e uma area que investiga o comportamento de materiais que admitem carac- 
teristicas elastica e de viscosidade quando submetidos a forgas de deformagao. No estudo destas 
caracteristicas, dois conceitos sao importantes: a tensao {stress) a{t) e a deformagao {strain) e{t), 
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ambas descritas como uma fungao do tempo t e os modelos matematicos aplicados procuram, 
justamente, descrever a relagao entre estas quantidades0 

Classicamente, as leis da mecanica que sao usa-das para descrever as relagoes supracitadas sao 
as de Newton e de Hooke, sendo a primeira para o comportamento de liquidos ideals (ou liquidos 
Newtonianos) modelado pela equagao 


o-W = (28) 

onde rj e a chamada constante de viscosidade do material. Ja a segunda lei (Hooke) modelada pela 
equagao 

a{t)=Ee{t), (29) 

descreve o comportamento de materials elasticos (ideals), sendo E a chamada constante elastica 
(ou modulo de elasticidade) do material. 

No estudo de sistemas modelados pelas equagoes Ea. ipSl) ou Ea. d^ . o interesse e obter des- 
crigoes sobre as respostas da tensao com res-peito a deformagao ou viceversa, sendo que os ensaios 
experimentais usualmente feitos sao os de [S]: 


• Relaxagao de tensao via deformagao controlada: onde se mede o estado de tensao 
causado no material de teste pela aplicagao de uma deformagao previamente definida. Neste 
caso, a resposta da tensao a aplicagao da deformagao e chamada modulo de relaxagdo G(f). 

• Teste de fluencia via tensao controlada: onde se aplicam ao material teste uma tensao 
previamente definida e medem-se as deformagoes causadas por essa tensao aplicada. Neste 
caso, a resposta da deformagao a aplicagao de uma tensao e chamada de fluencia 3{t). 


No estudo destes experimentos, as excitagoes impostas ao sistema (ou material teste) sao usual¬ 
mente modeladas pelo impulse (fungao ge-neralizada) 5-Dirac ou a fungao degrau (ou de Heaviside) 
H (t) definida por 


H{t) = 


0, t < 0, 
1 , <> 0 , 


(30) 


sendo que esta ultima e fisicamente mais realista que a primeira. Logo, neste trabalho, conside- 
raremos as respostas G{t) e J(t) quando submetidas a uma excitagao do tipo fungao degrau H{t) 
(Ver Figura 1 e Figura 2). 

Na natureza, entretanto, nao existem solidos ou fluidos ideals de modo que na pratica os ma¬ 
terials reals possuem propriedades situadas entre estes dois casos limites. Esquematicamente, o 
comportamento elastico (ideal) e descrito por uma mola, enquanto o comportamento viscoso (ideal) 
e descrito como um amortercedor (Figura 3) e os modelos aplicados na pratica sao uma combinagao 
destes dois elementos ideals. 

Entre os modelos mais simples, estao os de Maxwell e de Voigt, como descritos pela Figura 4: 

No modelo de Maxwell, a equagao que des-creve esta configuragao e dada por 


1 , , Id , , d , , 

-a(i) + --a(t) = -,(f) 


(31) 


^Em geral, a tensao e a deformagao sao campos tensoriais, assim para esta apresentagao pressupomos que os 
materials sao isotropicamente uniformes. 
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Figura 1: O Modulo de relaxagao para um solido ideal {E = I; trago continuo) e um fluido ideal 
(?7 = 1; trago pontilhado) gerados a partir de uma excitagao com a fungao de Heaviside. 



-10 12 3 

t 


Figura 2: O Modulo de fluencia para um solido ideal {E = 1; trago continuo) e um fluido ideal 
(?7 = 1; trago pontilhado) gerados a partir de uma excitagao com a fungao de Heaviside. 



(a) (I)) 


Figura 3: (a) Elemento elastico (ideal) e descrito por uma mola; (b) enquanto o elemento viscoso 
(ideal) e apresentado como um amortecedor. 
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Figura 4: (a) Esquematizagao em serie (Maxwell); (b) Esquematizagao em paralelo (Voigt). 


e as fungoes de resposta (modulo de relaxagao Gufit) e fluencia JM(t)) sao descritas como 

GM{t)=Ee-^\ (32) 

(33) 


Jjwii) — ~^E[(t) H— t. 

h Tj 


Ja no modelo de Voigt, a equagao que des-creve esta configuragao e dada por 


a(t) = Et(t) + i]—e(t) 


(34) 


e as fungoes de resposta (modulo de relaxagao Gy(t) e fluencia Jy(t)) sao descritas como 


Gy(t) = EH(t) + ??(5(t), 
^v(t) = ;^ (i - . 


(35) 

(36) 


Estes dois modelos tambem diferem das observagoes experimentais em alguns pontos de modo 
que estudiosos da area elaboram os mais diversos modelos (e.g., Zener, Kelvin, Burger, etc...) 
para tentar chegar a uma descrigao com maior sintonia com os experimentos [9]. Por exemplo, 
continuando com a tentativa de ge-neralizagao dos modelos de Maxwell e Voigt, os modelos de 
Zenei@ combinam os dois anteriores dos seguintes modos (Figura 5): (a) o modelo de Maxwell 
modificado combina um elemento de mola com o elemento de Maxwell em paralelo (b) e o modelo 
de Voigt modificado combina um elemento de mola em serie com o elemento de Voigt. 

O modelo de Maxwell modificado e descrito pela equagao 


T] di 




(37) 


^Seguindo a sugestao de nomenclatura como proposto em [^, estes modelos foram estudados por Clarence Melvin 
Zener (1905-1993) e os resultados publicados no sen livro de 1948, Elasticity and Anelasticity of Metals, Univ. of 
Chicago Press; e tambem sao conhecidos como modelos de solidos lineares padrdo. 
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Figura 5: Os dois tipos de modelos estudados por Zener: (a) Maxwell modificado; (b) Voigt 
modificado. 


com respostas (modulo de relaxagao GzAiit) e fluencia Jzm(O) dadas por, respectivamente 


^2 * 

GzAiit) = El + E2e ’> , 


E 1 E 2 

1 £'2^ V{Ei+E2) 

3zM{t) = {^Ei + E 2 ) ■ 


(38) 

(39) 


Ja o modelo de Voigt modificado e descrito pela equagao 


El + E 2 

V 


a{t) 


d(T(t) E1E2 


dt 


e{t) + El 


d£(0 

di 


(40) 


com respostas (modulo de relaxagao Gzv{t) e fiuencia 3zv{t)) dadas por, respectivamente 


Gzv{t) 

3zv{t) 


E? E1+E2 f E1E2 

-e H-, 

El + E2 El + E2 

El + E 2 1 -Ei-t 

-e ’> . 

El E 2 E 2 


(41) 

(42) 


A Figura 6, mostra as respostas dos quatro casos mencionados: Maxwell, Voigt, Maxwell mo¬ 
dificado e Voigt Modificado. Observamos que no modelo de Maxwell, o processo de relaxagao da 
tensao parece razoavelmente realista, no entanto, a resposta de fluencia mostra uma deformagao 
crescendo indefinidamente. No mo-delo de Voigt, o contrario acontece, a resposta de fiuencia parece 
razoavelmente realista, mas o processo de relaxagao e inexistente (permanece constante). Por outro 
lado, nos dois modelos modificados as respostas sao, ate certa forma, semelhantes, sendo que ambos 
qualitativamente ficam de acordo com alguns resultados experimentais, mas ambos possuem pro- 
blemas na des-crigao quantitativa, visto que predizem um valor finite para a tensao no tempo inicial 
t = 0 mesmo diante do salto (descontinuo) de mudanga de deformagao causada por uma excitagao 
do tipo Heaviside H{t) e nenhum dos dois atinge uma relaxagao de tensao completa. Alem disso, 
ambos modelos modificados predizem uma descontinuidade da resposta de fluencia imediatamente 
apos a aplicagao de uma excitagao de tensao do tipo Heaviside H(t). Estudos indicam que os mode¬ 
los podem ser aperfeigoados se forem incluidos mais elementos (de Newton, de Hooke, de Maxwell, 
de Voigt, etc) em serie ou paralelo, mas isto implicaria em equagoes cada vez mais complicadas 
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(!>) 


Figura 6: (a) Respostas do modulo de relaxagao da tensao para os modelos de Maxwell, Voigt, 
Maxwell modificado e Voigt modificado. (b) Resposta da fiuencia da deformagao para os 
respectivos modelos. 


e de ordens superiores. Uma alternativa para encontrarmos uma solugao mais satisfatoria a este 
problema, faz uso do calculo fracionario no equacionamento do problema como veremos a seguir. 

O leitor pode nao ter percebido, mas observamos um aspecto interessante a respeito dos dois 
casos ideais limites: a lei de Hooke idealizada pela Ea. (l2^ nos diz que a tensao a{t) e diretamente 
proporcional a derivada de ordem zero da deformagao e(t); enquanto a lei de Newton idealizada pela 
Eg. (1281) nos diz que a tensao a{t) e diretamente proporcional a derivada de ordem um da deformagao 
e{t). Entao uma generalizagao para o comportamento viscoelastico dos materials, segundo a ideia 
de derivadas de ordem fracionaria, seria conjecturar que a tensao a{t) e diretamente proporcional 
a derivada fracionaria de ordem a, sendo 0 < a < 1, segundo a equagao do tipc@ 

a(t) = Er“D“e(t), T=|, (43) 


sendo rj e E os coeficientes de viscosidade e elastico do material, respectivamente. Claramente 
os casos ideais limites sao recuperados quando tomamos a = 0 (Hooke) ou a = 1 (Newton). O 
operador seria qualquer operador de diferenciagao fracionaria escolhido para ser usado (e.g., 
Riemann-Liouville, Caputo, etc...). Em parti-cular, se quisermos um modelo que inclua todo o 
historico do material (em fungao do tempo), podemos estabelecer que o limite inferior seja t = —oo 
e ainda, exigindo a causalidade do sistema (i.e., que as fungoes sejam identicamente nulas para 
t < 0), entao as definigoes segundo Riemann-Liouville e Caputo coincidiriam. 

As fungoes de respostas (modulo de relaxagao GsB{t) e fiuencia Jsb(<)) relacionadas ao mo-delo 
da Ea. (|43)l podem ser obtidas, por exemplo, via a metodologia da transformada de Laplace como 
discutido na segao anterior e resultam em 


GsB{t) = 


Et°^ 




JsB(t) 


1 

Et°‘T (1 -f v) 


(44) 

(45) 


Estudos comprovam que mesmo este modelo simplista, resulta em descrigoes qualitativas mais 
precisas do comportamento de materials reals quando submetidos aos mesmos testes descritos na 
segao anterior e mesmo as descrigoes quantitativas parecem mais promissoras. De fato, podemos 


®Este modelo e conhecido como modelo de Scott-Blair [^. 
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(•) 


(b) 


Figura 7: Respostas de modulo de relaxagao para o modelo fracionario de Scott-Blair para 
diferentes parametros a, {E = 1, r] = 1); (b) Respostas da fluencia da deformagao para o modelo 
fracionario de Scott-Blair para diferentes parametros a {E = 1, rj — 1). 


observar pela Figura 7, que todas as respostas (modulo de relaxagao GsB{t) e fluencia Js_B(t)) para 
diversas ordens a se comportam quantitativamente como o esperado, i.e., o modulo de relaxagao 
possui um valor infinitamente grande para t = 0 e decaindo para relaxamento total quando t —>■ oo, 
ja a fluencia mostra o crescimento da deformagao de forma mais apropriada com os experimentos 
mantendo a continuidade em t = 0, lembrando que os graficos dados na Figura 7 foram feitos a 
partir de excitagoes modeladas pela fungao degrau 77(t). 

Este exemplo, mesmo que simples, mostra o potencial do uso do calculo fracionario para abor- 
darmos o problema de viscoelasticidade linear e de fato, este ramo tern se mostrado um terreno 
fertil para estudiosos da area com diversos livros [ziiniiii] mostrando detalhamento dos modelos e 
resultados experimentais. 

Mencionamos tambem outras areas em que o calculo fracionario tern se mostrado uma ferramenta 
poderosa, a saber, na descrigao de osciladores harmonicos fracionarios, forgas de fricgao fracionarias, 
controladores fracionarios, equagoes de ondas fracionarias, etc... [iiniiiiiMiiiT]. 

5.2 Osciladores Harmonicos 

Nesta segao discutimos uma possfvel genera-lizagao do problema associado ao oscilador harmonico. 
Recuperamos o caso inteiro e a esse propomos uma generalizagao no sentido fracionario, isto e, 
introduzimos dois parametros associados a ordem da equagao fracionaria. A fim de obtermos 
uma solugao dessa equagao, utilizamos a metodologia da transformada de Laplace e fornecemos 
essa solugao explfcita para o caso em que nao temos o termo associado com o atrito, ou seja, 
consideramos o oscilador sem amortecimento e apresentamos essa solugao em termos das fungoes 
de Mittag-Leffler. 

Vamos considerar apenas o caso em que a derivada e tomada no sentido de Caputo (Definigao 
n enquanto as condigoes iniciais admitem a classica interpretagao. Para discutirmos tal problema 
inicia-se por obter uma equagao integral correspondente a equagao diferencial associada ao oscilador 
harmonico 

d2 d 

^x(t) + + uj^x(t) = 0 

com /J. > 0, representando o termo associado ao atrito e w > 0 e a frequencia do oscilador, sendo as 
condigoes iniciais a;(0) e x'(0) dadas. 
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Integrando essa equagao diferencial duas vezes podemos escrever 

x(t) = x{0) + /itx(O) + tx'{0) — 

pt pt pV 

^ / x(t;)du—/ / x(M)du dx 

Jo Jo Jo 

que, apos a utilizagao do teorema de Goursat laiio] permite escrever 

x{t) = x(0) + x(0) + t x'(0) — 

/i / x{y)Av — (jj^ / {t — u)x(u)Au. 

Jo Jo 

Entao, dados x(0) e x'(0) obtemos uma equagao integral associada ao problema do oscilador 
harmonico equivalente ao problema composto pela equagao diferencial e as condigoes iniciais. 

A fim de explicitar os calculos, consideremos o oscilador harmonico livre, isto e, sem atrito com 
a equagao integral ja na forma fracionaria 


xit) = x(0) + tx'(O) — / (t — drt 

r(a) Jo 


(46) 


com 1 < a < 2. Escolhemos esse intervalo a fim de recuperar o resultado do oscilador harmonico 
classico, no caso em que a = 2. 

Note que essa equagao integral fracionaria corresponde a equagao diferencial fracionaria 

f + a;“x(t) = 0 

e as condigoes iniciais x(0) e x'(0) sendo 1 < a < 2 e a derivada considerada no sentido de Caputo. 

(Definigao [3]). 

Vamos procurar uma solugao da equagao integral atraves da metodologia da transformada de 
Laplace que apos aplicada na Ea. (l46t e utilizando a definigao do produto de convolugao permite 
escrever 


s Jo r(a) 


-kx{t)e 


onde o simbolo * denota o produto de convolugao. Logo, obtemos 

s 


sendo 


poo 

F{s) = / x{t)i 
Jo 


dt 


a transformada de Laplace de x{t) de parametro s com Re(s) > 0. Isolando F(s) podemos escre-ver, 
ja rearranjando, 


F{s) = x(0)- 


r'(0): 


Para recuperar a solugao x{t) tomamos a transformada de Laplace inversa. Logo para 

nC-\-i00 




e®‘F(s)ds = £-i[E(s)], 
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obtemos a solugao na forma 

x{t) = x(0) Ec{-uj°‘r) + a;'(0) Sa, 2 (-w“f“) 

onde Ea{-) e Ea^p{-) sao as fungoes de Mittag-Leffler com um e dois parametros, respectivamente 
liiniin]. Admita, enfim, que a:(0) = 1 e x'(0) = 0 logo a solugao neste caso e 

x{t) = Eci{—uj°'t°‘) 

que, no caso extreme a = 2 fornece x{t) = cos cot que e a solugao do problema associado ao oscilador 
harmonico de ordem inteira. 


6 Conclusao 

Neste trabalho foram introduzidas as chamadas equagoes diferenciais fracionarias e constatado que 
a metodologia da transformada de Laplace mostra-se, assim como no caso das EDOs de ordem 
inteiras com coeficientes constantes, uma otima ferramenta a ser usada em vista dos benefi'cios da 
sua simplicidade. Em particular, verificou-se de forma indireta que da mesma forma que a fungao 
exponencial desempenha um papel de destaque no calculo classico de ordem inteira, as fungoes 
de Mittag-Leffler desempenham um papel destacado analogo no calculo de ordem arbitraria. O que 
queremos denotar por esta afirmagao e que assim como se faz necessario um bom entendimento da 
fungao exponencial para o estudo e aplicagoes do calculo diferencial e integral classicos o mesmo 
pode ser dito a respeito das fungoes de Mittag-Leffler para aqueles que desejam se aventurar pelo 
calculo fracionario. De fato, nao e incomum ao resolvermos equagoes diferenciais (ordinarias) de 
ordem fracionaria, obtermos solugoes que, se nao diretamente expressas em termos de uma fungao 
de Mittag-Leffler, entao ao menos se relacionam com ela atraves de alguma identidade. 

Por fim, observamos que, de modo relativamente simples, problemas classicos da Fi'sica-Matematica 
ficam bem modelados sob a optica do calculo fracionario e dao resultados promissores (inclusive 
coincidindo com os resultados de ordem inteira ja firmados) mesmo quando sao elaborados de forma 
simples, como verificados nos exemplos de viscoelasticidade linear e osciladores harmonicos. 

O leitor interessado em conhecer mais modelagens usando o calculo fracionario, pode recorrer 
as fontes [l2l[Mll29l[30l[3ll|3l[33l|3l|35]. 
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